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ABSTRACT 



The dust size distribution in molecular clouds can be strongly affected by ice-mantle formation and (subsequent) grain coagulation. 
Following previous work where the dust size distribution has been calculated from a state-of-the art collision model for dust aggregates 
that involves both coagulation and fragmentation (Paper I), the corresponding opacities are presented in this study. The opacities are 
calculated by applying the effective medium theory assuming that the dust aggregates are a mix of 0. 1 yum silicate and graphite grains 
and vacuum. In particular, we explore how the coagulation affects the near-IR opacities and the opacity in the 9.7 fim silicate feature. 
We find that as dust aggregates grow to yum-sizes both the near-IR color excess and the opacity in the 9.7 yum feature increases. Despite 
their coagulation, porous aggregates help to prolong the presence of the 9.7 jim feature. We find that the ratio between the opacity 
in the silicate feature and the near-IR color excess becomes lower with respect to the ISM, in accordance with many observations of 
dark clouds. However, this trend is primarily a result of ice mantle formation and the mixed material composition of the aggregates, 
rather than being driven by coagulation. With stronger growth, when most of the dust mass resides in particles of size ~10 pm or 
larger, both the near-IR color excess and the 9.7 jjm silicate feature significantly diminish. Observations at additional wavelengths, in 
particular in the sub-mm range, are essential to provide quantitative constraints on the dust size distribution within dense cores. Our 
results indicate that the sub-mm index /3 will increase appreciably, if aggregates grow to ~100 /jm in size. 

Key words. ISM: dust, extinction - ISM: clouds - Opacity - Infrared: ISM - Submillimeter: ISM 



1. Introduction 

Interstellar grains are an important component of molecular 
clouds. Interstellar dust is a major source of opacity in the 
far-ultraviolet region of the spectrum and hence the dust char- 
acteristics are an important aspect of the spectral energy dis- 
tribution of molecular clouds. In addition, through their in- 
fluence on ph otodestruction rates, dust controls the survival 
of molecules dRoberge et al.l Il99ll) . Grain surfaces also pro- 
vide a 'meeting' place for atomic and mole cular species and 
hence enable reactions to pr oceed efficiently dTielens & H agen 
1982; Hasega wa et al.l [T992). Thus, molecular hydrogen is gen- 
erally thought to be formed 'exclusively' on grain surface 
dHollenbach & Salpeterll97ll) and the formation of other species 
such as w ater may also well be dominated by gra in surface 
chemistry (ICharnlev et al.|[l992t lHasegawa et al.ll 1 992b . Finally, 
grains are often used as a proxy for the gas and observations at, 
in particular, sub-millimeter wavelengths pr ovide a convenient 
probe of the density distribution of clouds dJohnstone & Bally 
I2006T) . Such observations are therefore often employed in ana- 
lyzing the star formation characteristics of molecular clouds in 
terms of, for exam ple, the masses of cloud cores and prot ostel- 
lar condensations (lAlves et alj2007t lJ0rgensen et all "2008). The 
physical characteristics of interstellar dust are therefore key to 
understanding many aspects of molecular clouds and their evo- 
lution. 

The properties of inters t ellar d ust, however, evolve in - 
side dense clouds dDwekl Il998t IZhukovska etall 120081) . 



Observationally, large grains are indicated inside dense clouds 
through an increase in the value of the total-to-selective extinc- 
tion ratio, Ry, from about 3.1 in the diffuse interstellar medium 
to values as high as 5.5 for sight lines traversing dense core s 
(IWilking et alJl980t IWhittetH2005T:IOlofsson & 01ofssonl20Toh . 
This is supported by studies of the wavelength dependence of 
polar ization which also indicates grain growth dCarrasco et al.l 
1973). The unusually low value of the visual extinction per 
H-atom towards the dense cloud, p-Oph implies that grain 
growth - at least for this cloud - reflects coagulation rather 
than growth through mantle formation (Jura 1980). Early stud- 
ies have revealed that for low column densities in the diffuse 
ISM, the near- and mid-IR extinction shows a 'universal' behav- 
ior characterized by a power-law dep endence of extinction on 
wavelength with an exp onent of -1.8 (Rie ke & Lebo fskv 1985; 
Martin & Whittet 1990). Recently, large scale surveys at near- 
and mid-IR wavelength (2MASS, UKIDSS and GLIMPSE) have 
provided new insights into the extinction behavior of dust deep 
inside dense clouds by probing the colors of background stars. 
These studies have revealed that the near-IR extinction law 
shows a flattening in dense cores, ty pically at depths corre- 
sponding to an Ay of ~20 magnitudes (jlndebetouw et al.l f2005; 
IChapman et al. 2009; Cambresv e t al.ll2.Ql H) indicative of grain 
growth dWeingartner & D raine 2001). Recently, too, near-IR 
scattered light has been detected from dense cloud cores and this 
is indicative of grain g rowth to micron-sizes (Pa gani et al 1 2010: 
Steinacke r et al.ll2010l) . Finally, grain growth in dense clouds is 
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also supported by an analysis of the sub-millimeter emi ssion 
( Bia ncni etal.ll200l iKramer et al.ll2003b LShirlev et al.ll2011h . 

The high sensitivity of the Spitzer/IRS instrument has al- 
lowed for the first time studies of the profile of the 10 pm sili- 
cate feature for dust inside dense clouds by probing background 
stars. These studies have revealed two new aspects of the ex- 
tinction behavior of dust inside dense clouds. First, while the 
ratio of the mid-IR extinction to the near-IR extinction increases 
with increasing depth (the flattening of the extinction curve men- 
tioned above), the total extinction at 10 yum does not increase; 
e.g., per unit near-IR extinction, the increase in the 10 pm 'con- 
tinuum' extinction is accompanied by a concomitant decrease in 
the strength of the silica te feature (IChiar et alJ 120071: iMcClurd 
l2009t [Boogert etaDl201lh . This decrease in the optical depth of 
the silicate band relative to the total near-IR extinction is accom- 
panied by a s hift to the blue in the profile compared to th at of the 
diffuse ISM (van Breem en et al ]l201lHChiareFaT]l2011h . In this 
study we will investigate the role of dust coagulation on the be- 
havior of the 10 /mi silicate opacity vs. the near-IR color excess. 

Variations in the dust properties inside dense clouds can re- 
flect accretion of ice mantles and/or coagulation. Taking the 
measured strength of the absorption features due to interstel- 
lar ice mantles towar ds background sources behind dense clouds 
(Boo gert et al.ll20lTl) and the intrinsic strength of these bands, 
the increase in grain volume due to ice is ^ 2.5 x 10~ 27 cm 3 per 
H-atom (cf. iTielens & Allama ndola 1987); a little less than the 
silicate dust volume (^3.6 x 10 -27 cm 3 per H-atom). However, 
these ice mantles only increase the extinction within specific ab- 
sorpti on modes and have little infl uence on the overall extinc- 
tion (Ossenko pf & Henning| [r994). Moreover, because the sur- 
face area of grains is dominated by the small end of the size dis- 
tribution, the mantle thickn ess is very lim ited (< 175 A if all oxy- 
gen were to deplete as ice ; iDrainell 19841) . In terms of grain size, 
c o agul ation is the pore important one of the grain growth pro - 
cesses dOssenkopjl 1 9931: IWeidensc hilling & Ruzmai kinall 1 9941) . 
At low collision velocities, small dust grains will stick upon 
collision forming larger conglomerates. Collisions a mong these 
aggre gates builds then larger and larger grains (cf. Ossenkopf 
1993). Because the coagulation process itself as also the ex- 
tinction properties depend strongly on the grain structure, much 
theoretical and experimental effort has focused on understand- 
ing the evolution of the porosity of co l hsional aggregates 
(Ossenkopf 1993; Dominik & Tielens 1997; Paszun & Dominik 
2009: iBlum & Wurml 120081; ISuvama et alJ 120081; IWada et al.1 
2008; lOkuzumil l2009h . These studies show that initially coag- 
ulation may lead to the formation of open, fractal structures, but 
eventually, if collision velocities b ecome high, comp action and 
then fragmentation will take over (lOrmel et al.1120091) . As a re- 
sult, the structure of the aggregates - and therefore the extinc- 
tion properties - depend on the density and history of the cloud. 
Recently, we have calcula ted the colhsional growth of interstel- 
lar grains in dense clouds (lOrmel et alJ2 009). Here, we examine 
the implications of this growth for the extinction behavior of dust 
in dense clouds. 

This paper is organized as follows. Section|2] provides an 
overview of previous work and details the assumptions em- 
ployed in the calculation of the opacities. In Sect. [3] the extinc- 
tion coefficients are calculate for solid grains, varying only their 
composition (silicates, graphites, ice-layers, and mixed vari- 
ants). Sect. [3] also discusses the sensitivity of the porosity of 
the dust aggregates. In Sect. [4] we apply the opacity calcula- 
tions to the aggreg ate dust size distribution obtained from the 
lOrmel et"aTI ((2009) study. Section[5] discusses several implica- 
tions and summarizes. 



2. Overview of previous work and methodology 

The dust opacities that we present in this work are the result of a 
chain of modeling efforts. These include: 

1 . a colhsional dynamics code to perform a parameter study on 
the outcome of a colhsion between two aggregates of grains 
dPaszun & Dominikll2009l) : 

2. a colhsional evolution code to compute t he size distribution 
of dust aggregates as function of time (lOrmel et alJ 2009; 
henceforth Paper I); 

3. an effective medium approach t o calculate the op tical prop- 
erties of porous dust aggregates (Min et al. 2008); and 

4. a Mie scattering code. 

In the following, we briefly review these steps. 

2.1. Collisions among dust aggregates 

At low velocity, collisions among dust grains result in loosely 
bound structures held together by surface forces that act on the 
contact points between the individual grains. These structures 
are referred to as dust aggregates and subsequent 'growth' in- 
volves collisions between two dust aggregates. The outcome of 
a collision between two dust aggregates depends on the material 
properties, the internal structure of the aggregates, and on the 
collision properties (impact parameter, collision velo city, mass 
ratio). In a pioneering work, Dom inik & Tielensl d 1997b derived 
collision recipes to predict the outcome of collisions between 
two arbitrary aggregates, using appropriate scaling behavior. In 
particular, the outcome of a collision between two grains or ag- 
gregates - sticking, compaction, or fragmentation - depends on 
the collision energy in comparison to a critical energy and the 
numb er of contacts within the aggregate (see IBlum & Wurml 
120081 for a review). 

Using a 3D molecular dynamics code, in which the equa- 
tions of m otion for each indiv i dual g rain within the aggregate 
ar e solved. iPaszun & Dominikl J2009) have expanded the study 
of Dominik & Tielens (1997) by including fragmentation and 
the influence o f porosity on the collisio n outcome. Using look- 
up tables, the IPaszun & Dom inik (2009) recipe provides a pre- 
scription for the outcome of the collision: the size distribution 
of the colhsional fragments (when present) and of the porosity 
increase/decrease of the collision products. The great advantage 
(for modelers) of using the recipe approach is that we only need 
to follow the average properties of the particles (the porosity and 
size of the aggregates, material properties of the grains) and the 
collision properties (i.e., the velocity field), and can dispense of 
modeling the full substructure of dust aggregates. Indeed, the 
latter approach is computationally too expensive since molecu- 
lar dynamics code can only treat a limited number of grains. 

2.2. Application to molecular cloud cores 

In Paper I we have included these recipes in a Monte Carlo 
code to compute the size distribution of dust aggregates as func- 
tion of time. In Monte C arlo codes (e.g., lOrmel et"aL1 12007; 
IZsom & Dullemondll2008l) the computational particles are char- 
acterized by their properties. In this case, the aggregates were 
represented by their mass m and their projected surface-area over 
mass ratio, cr/ra. An accurate prescription for the latter property 
is very important, as it plays a critical role in the determination 
of the relative velocities among the dust aggregates. In Paper I 
we explicitly follow both <x and m, i.e., cr is not a priori assumed 
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Fig. 1. The size distribution of dust aggregates as function of time as calculated bv lOrmel et alj d2009l) for bare grains (left panel) 
and ice-coated grains (right panel). The gas density is n g = 10 5 cm- 3 . For bare grains the size distribution evolves towards a steady 
state due to the emergence of fragmentation, while aggregates with ice-coated grains keep growing. The gray shading denotes the 
spread obtained from independent runs of the Monte Carlo code. The lower panels show the filling factor of the dust aggregates. 
Decreasing filling factors indicate a fractal structure. But for larger aggregates compaction halts the fractal growth. 



to be a function of m as, e.g ., for compact sphe res (cr <x m ' 3 ) or 
very fluffy particles cr oc m (Minato et al. 200^). 

Equivalently, we can express cr as a size, a a = ycrfn, and 
define the geometrical filling factor as 

Volume occupied by N grains _ N x 47ra 3 /3 _ /Va 3 
°" Volume equivalent to surface area 47ra 3 r /3 aj T 

(1) 

where we have assumed that all N grains are spherical and of the 
same size ao. This monodisp erse grain approximation i s used in 
both the parameter study of Paszun & Dominik ( 2009) as well 
as in Paper I; in the latter study the radius of the grains ao was 
fixed at 0. 1 //m. 

In this study we furthermore limit our results to models that 
start out at a gas density of n g = 10 5 cirT 3 . Figure Q] shows the 
evolution of size distribution for this density. Here, f(N) gives 
the number density spectrum of aggregates by their number of 
grains, N. Following Paper I, we multiply by A^ 2 to obtain a 
proper comparison for the mass density for a logarithmic scaling. 
We expect that for other gas densities the coagulation will fol- 
low a similar trend (at least initially) but at different timescales 
(see Paper I). Paper I considered two modes of aggregation: ag- 
gregation among bare (silicate) grains and among ice-coated 
grains. For the sticking behavior, it is the surface properties that 
matter and even a modest amount of freeze out on bare (sili- 
cate or carbonaceous) grains can significantly speed-up the co- 
agulati on, since the critical energies are expected to be much 
larger (Dom inik & Tielenslll997l: IWada et alj r2009). We indeed 
observed a rather sensitive dependence of our results on these 
material properties. In the bare silicate models the aggregates 
were weak and were prone to fragmentation. Consequently, a 
steady state between coagulation and fragmentation with little 



net growth was quickly established on timescales larger than 
~3 x 10 5 yr (Fig.QJ). However, if the grains were coated by a 
(thin) ice layer, small grains were quickly swept into aggregates. 
These aggregates continued to grow through mutual collisions. 
Initially, the aggregates remained rather porous as the collisional 
energies were low, but with increasing size and increasing en- 
ergy the aggregates started to compact. Only on timescales >10 7 
yr were collisional energies powerful enough to provide some 
form of erosion and replenish the smaller grains (Fig.QJi). 

In Fig.Q]:,d the (mean) filling factor (Eq. ([TJ) of the aggre- 
gates is plotted, as function of their size. Most of the curves fall 
on top of each other, indicating that (f> lT is a function of size only. 
Initially, the filling factor decreases with size, which indicates 
that the growth is fractal. Due to subsonic turbulent motions ag- 
gregates move at an appreciable relative velocity of «10 m s 1 
which increases with the growth of these particles. As a result, 
energetic collisions will halt the decrease in 4> a through com- 
paction. After t = 10 7 yr there is a noticeable increase in <p lT 
towards <p lT = 33%, whi ch is the maximum adopted value of <p lT 
for evolved aggregates (Guttler et al. 2009, Paper I). 

2.3. Optical properties of dust aggregates 

From the size distributions of aggregates provided by Paper I we 
compu t e opacities , follow ing the method described in lMin et"al1 
(2008l). lMin et alJ d2008) describe an efficient method to obtain 
the optical properties of porous dust aggregates using an effec- 
tive medium approach. This involves applying a mixing rule that 
appropriately averages the optical constants of the various com- 
ponents out of which the aggregate is constituted. The resulting 
effective refractive index (m e g) can then be u sed in Mie's solu- 
tion of light scattering. The iMin et ail d2008h study includes a 
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Fig. 2. The near/mid-IR extinction for spherical grains. The 
curves show the opacities corresponding to 0.1 yum spheres of 
silicate (light gray) and graphite (black). Dark gray curves corre- 
spond to a 1 :2 mix of graphite to silicates of 0. 1 //m grains (solid 
gray curve) and for the MRN distribution (dashed curve). The 
dotted curve shows the adopted continuum for calculating k s u, 
the amount of extinction in the 9.7 fim silicate feature. Shown 
for illustrative purposes is the observed ISM extinction curve 
(thin solid curve) with arbitrary scaling in the ^-direction. 

prescription to add vacuum as a separate component to the mix- 
ing rule. Here we app l y the m ethod to spherical grains. We con- 
sider the iBruggemanl (119351) mixing rule, with N c components 
and vacuum: 



fm^fiCHcOnt/mes) + (1 - fm)as c (l/m e s) = 0, 



(2) 



;=1 



where m e g is the effective refractive index, fg\\ the filling factor 
of the aggregates (see below), fi the volume fraction of com- 
ponent i with c fi = 1, m, the refractive index of component 
i, and a c (m) the polarizability of the components. For homoge- 
neous spheres, the polarizability is 



ahs(m) = 4na~ 



1 



m 2 + 2' 



(3) 



where a is the aggrega te's radius. For ice -coated spheres, the 
polarizability becomes (Ivan de Hulsdfl98lh 

gcs(g) ^w (£2 - i)(e ' +2g2)+e ; (2g2+i)(fii - fi2) , ( 4) 

(e 2 + 2)(ci + 2e 2 ) + Q\2s 2 - 2)(e, - e 2 ) 

where si = m 2 is the dielectric constant of the inner core, £2 
the dielectric constant of the coated material, and 1 - Q the frac- 
tional size of the coated material; that is, the inner core consists 
of fractional volume g 3 - In our case, £1 represents silicates or 
graphite and £ 2 ices, where we take Q _ 0.9 for a t hin coating. 
The refractive indices are taking from fDraind ( 12003b . 

Thus, Eq. (f2]l is solved for m e g (or £ e ff). In general, this re- 
quires us to solve a polynomial equation, which can be done 
numerically or analytically if the degree is lo w. The only other 
remaining parameter is the filling factor fm . iMin et a I1 (l2008h 
suggests to use the radius of gyration of the aggregate a g , i.e., 
/en = Na^/ag, and reports an excellent match with more rigor- 
ous coupled dipole approximation calculations. However, from 
Paper I only a ff (rather than a g ) is available, since we have not 



computed the full substructure of the aggregates. For this study 
we will assume that these radii are similar, a a w a g . For suffi- 
ciently dense aggregates of fractal dimension Df > 2, one can 
expect that the projected surface area is indeed similar to, al- 
though always smaller than, the outer-most radius (see Paper I 
for a quantitative relation between these radii). The same argu- 
ments holds for the gyration radius a g . For fluffy aggregates of 
Df < 2 these assumptions will fail, however, and our method 
becomes invalid. But in our case aggregates never become very 
fluffy (<f> always exceeds 10%, see Fig.[TJ. Therefore, these radii 
are similar, which allows us to substitute (f> lT for fm- 

From the effective optical constant m e g = n e s+ikes and a par- 
ticle of effective size a^, we apply the Mie scattering solution to 
Maxwell's equations (e.g., iBohren & Huffman| [l983). We then 
obtain the scattering cross section C sca (A), the absorption cross 
section C- d b s , and the extinction cross section, C ext = C sca + C a b s , 
as function of wavelength A. From these cross sections and 
the particle distribution function we finally obtain the mass- 
extinction coefficient, K ext (A), which gives the cross section for 
extinction per unit mass dust. 



3. Opacities for static grain distributions 

3.1. The near-IR extinction law 

TableQ] gives an overview of all opacity calculations. We first 
consider opacities resulting from individual grains, or grain mix- 
tures. In Fig.|2]the light gray curve gives the mass-weighed ex- 
tinction coefficient for 0. 1 fim silicate grains, the black curve cor- 
respond to graphite grains, and the dark gray curve corresponds 
to a silicate-graphite mix of grains in a 2:1 ratio (we motivate 
this choice for the ratio below in Sect. 13.2b . The 9.7 fim feature 
is prominently visible in the silicate curve and in the silicate- 
graphite mix: the silicates are the carrier for the 9.7 /mi feature. 
At near-IR wavelengths the carbonaceous grains dominate the 
opacity; they are the carriers for the E(J - K s ) extinction. 

For comparison, we also plot in Fig.|2]the observed ISM ex- 
tin ction curve (th i n solid line). At short wavelength s we adopt 
the | Cardelli et alj {1 989) fit (with coefficients from Fitzpatrick 
1999 and Ry - 3.1), whereas a t wave lengths >2.2 fim we adopt 
the profile of IChiar & Tielensl d2006h . The fit plotted in Fig.H 
provides A{A) multiplied by an arbitrary constant (only the shape 
is meaningful). The comparison with the observed ISM extinc- 
tion curve is meant to be for illustrative purposes; it is not our 
intention to actually fit this curve. 

W e have also calculated opacities for an MRN dMathis et al.l 
[1971 distribution of grains, where the grain size a is distributed 
according to a -3.5 power-law between a lower size of 50 A and 
an upper size of 0.25 /mi. The extinction coefficients are again 
calculated for a 2:1 silicate:graphite mix (dashed gray curve in 
Fig- [3- Comparing K ext to the 0.1 yum mix (solid gray curve) 
one recognizes minor differences. This is understandable since 
in the Rayleigh limit A- ext depends only on the mass of the grains. 
However, for A < 0. 1 /urn the opacities of the MRN distribution 
become larger than those of the 0. 1 fim grains due to their larger 
total surface area per unit mass. 

In Fig. [3] we illustrate the opacities resulting from the 
Bruggeman mixing rule Eq. (f2]i for a 3 [im silicate aggregate 
of different porosity, indicated by the filling factor fm- The fill- 
ing factor of the aggregates reaches a minimum of w 10% (see 
Fig.Q]). As explained above these are porous particles, but in no 
way resemble very fluffy aggregates of low fractal dimension. 
However, Fig.[3]shows that the effect on the opacities is not neg- 
ligible, especially at small wavelengths (A <sz 1 jum). For a filling 
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Table 1. Model runs. 



Name 


Mixing type 


Size distribution 


Components 


Figure ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


0.1 sil 




Spheres, 0.1 fim 


silicate 


Figs. |2] and |4] 


O.lgra 




Spheres, 0.1 fim 


graphite 


Figs. |2j and |4j 


O.lsil+gra 


I 


Spheres, 0.1 fim 


silicate, graphite 


Figs. |2| and |4| 


MRNsil+gra 


I 


Spheres, MRN 


silicate,graphite 


Figs. [2] and [4] 


U.llC-Sli 




Spheres, 0.1 


ic-silicate 


Flg-|4| 


O.lic-gra 




Spheres, 0.1 fim 


ic-graphite 


Fig.S 


O.lics+icg 


I 


Spheres, 0.1 fim 


ic-silicate,ic-graphite 


Fig. a 


0.1(sil,gra) 


II 


Spheres, 0.1 fim 


silicate,graphite 


Fig.S 


(sil,gra) 


II 


Aggregates 


silicate, graphite 


Fig.l6ji 


(ic-sil,ic-gra) 


II 


Aggregates, coated 


ic-silicate, ic-graphite 


Fig.[6j) 


(ic-sil,gra) 


II 


Aggregates, coated 


ic-silicate, graphite 


Fig.[6t 


ic-sil+gra 


I 


Aggregates 


graphite,ic-silicate 


Fig.p 



Note. — Overview of all model runs. (1) Model name. The top part of the table corresponds to the static gr ain distribut i ons co nducted in Sect.[3] 
whereas the bottom rows concern the opacity calculation of the (evolving) aggregate distributions from Orm el et all d2009h (Sect.S- (2) The 
considered mixing type (if applicable), see Fig. [5] Type I indicates the two components (Col. 4) are mixed in a spatial volume; type II indicates 
they are mixed within aggregates or grains. (3) Adopted size distribution. For the static grain distributions these are either monodisperse 0.1 
fim grains or follow an MRN distribution. For the aggregate models these corresponds to either Fig.QJ (non-coated) or Fig.fTJj (ice-coated), 
respectively. (4) Components included in the aggregate calculations. Here, 'ic' stands for ice-coated. The vacuum component that is present in the 
aggregate models is not listed. In every model where both silicates and graphites are present their mass ratio is fixed at 2:1. 




Fig. 3. Opacity changes due to porosity variation (indicated by 
the filling factor /gn) for a silicate aggregate of size 3 fim. 

factor of 10% the extinction is increased by a factor of 10. This 
effect can be understood from geometric considerations: the sur- 
face area per unit mass is increased by a factor 10. Moreover, 
the ~10 fim spectral feature is more pronounced than in the low 
porosity case. Porous aggregates preserve to a large extent th e 
spectral signature of their constituent grains dMin et al.l 12008). 
This also holds for the near-IR wavelengths; the fsn = 0.1 ag- 
gregates are optically more pristine. 

3.2. The silicate 9.7 fim and E{J - K) indicators 

We will now introduce two indicators that are used to assess 
the evolutionary state of a core. These are the strength in the 
9.7 fim silicate feature, obtained after the continuum subtraction 
(see Fig.[2| and the near-IR color excess, E(J - K) - Aj — Ak s , 
with A j and Ag s the extinctions in terms of magnitudes in the J 
(1.25 fim) and the K s (2.2 fim) bands, respectively. These quan- 
tities are plotted against each other in Fig. [4] where we applied 
the conversion between extinction in terms of optical depth or 



Fig. 4. The opacity in the 9.7 fim silicate absorption feature vs. 
near-IR color excess for spherical grains. The colors represent 
the material: open symbols represent silicates, black symbols 
represent graphite, and gray models represent a mix of 2/3 sili- 
cate and 1/3 graphite grains. Models involving ice-coated grains 
are indicated by circles and the diamond indicate an MRN dis- 
tribution of spherical grains. The 6) . 1 (sil , gra) model denotes 
the opacities of a 0. 1 fim sphere that consist for 2/3 out of sili- 
cates and 1/3 of graphite. The solid line of slope q = 0.34 corre- 
sponds to the diffuse ISM. is the total dust column in units of 
g cirT 2 . 



magnitudes, A, = 1.078r,. In Fig.|4]the extinction is normalized 
by the total dust column, (units: g cm -2 ); to obtain the total 
extinction (t s h or E(J - K)) for a particular core one must mul- 
tiply by its Y.J value. Clearly, for the pure silicate model (open 
square) the near-IR color excess is negligible, whereas for the 
pure graphite model (black square) the optical depth in the 9.7 
feature is unimportant. For diffuse ISM clouds, observations of 
the silicate optical depth and the near-IR color excess all collapse 
on a straight line with constant slope, q = T s g/E(J - Ks ) ~ 0.34 
(iRoche & Aitkenll984tlWhittell2003HChiar et al.l2007l) : the dif- 
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Fig. 5. The two mixing types considered in this study. In type I 
mixing aggregates consisting out of grains of the same material 
are mixed; the opacity is calculated for each of the components 
individually and averaged. In type II mixing the materials are 
mixed on the level of the individual aggregates, which adds an- 
other component to the Bruggeman mixing rule of Eq. (f2]). 



fuse ISM curve. The slope therefore measures the underlying 
opacity ratio between the indicators - a function of the grain size 
distribution and composition - which may be assumed constant 
in case of the diffuse ISM. 

To reproduce the diffuse ISM curve the silicate and car- 
bonaceous materials have to be mixed. We find that with a 2:1 
mix of silicates to graphite the opacity ratio falls approximately 
on the diffuse ISM curve (gray symbols). Therefore, we have 
adopted this ratio in all our calculations. In Fig. [4] we also plot 
the indicator-values for ice-coated grains (circles) for a coating 
of 10% in radius. This has the effect of reducing the strength 
of the E(J - K s ) colors by «25%. Therefore, an ice-coating can 
causes the indicators to deviate from the ISM value. 

Finally, we have calculated the opacities for a spherical grain 
that consist for 2/3 of silicate and 1/3 of graphite; i.e., we calcu- 
late the emergent opacity by applying Eq. (f2]i with N c = 2 and 
/(in = 1. This mixing of the material components on the level 
of the individual aggregate or grain is referred to as 'type II 
mixing' (see Fig. [5j and is indicated by the brace notation in 
Table[T] e.g., . 1 (sil , gra). In contrast, the mixing of the in- 
dividual components in space that we previously encountered 
is referred to as type I mixing and is denoted with a '+' sign, 
e.g., 0. lsil+gra. Figure[5]provides a graphical illustration of 
the two mixing types, in case of aggregates. From Fig. [4] it is 
seen that for type II mixing the emergent opacities are lower 
than for type I, mixing. In type I mixing the optical properties 
of the participants are conserved individually and the mixed op- 
tical properties are a linear average of the individual opacities. 
However, in the effective medium approach used with type II 
mixing (Eq. (f2]l) no such linearity is present! In addition, the 
continuum subtraction applied to determine t s ;i is also different 
between both cases. From these results, therefore, we anticipate 
that the formation of aggregates composed of different materials 
and (to a minor extent) the formation of ice-mantles causes the 
indicators to deviate from the diffuse ISM curve. 

4. Opacities for evolving dust aggregates 

4.1. Time-dependent extinction law 

The bottom rows of Table Q] list the four aggregate models for 
which we have performed opacity calculations. The correspond- 
ing IR-spectra are given in Fig. [6] For all these models, silicates 
and graphites are assumed to be present at a 2:1 ratio. In Fig. [6] 
the opacity distributions are shown for the same times as the 



size distribution in Fig.Q] We consider four aggregate models, 
where we vary the size distribution (Fig.QJ or b), and the com- 
position of the aggregates (material components, coating). In the 
first three models we assume that the aggregates consist of a 
mix of three materials: silicate, graphite and vacuum; i.e., N c in 
Eq. (f2]i equals 2 (type II mixing). However, in the last model we 
assume that the graphite and silicates are (spatially) separated 
and constitute different distributions. Then, we solve for each of 
the k v individually and average to obtain the mixed value (type I 
mixing). 

In Fig. [6^ the opacities corresponding to models without ice- 
coating, i.e., as shown in Fig.QJ, are presented (sil , gra). It is 
clear that the opacities evolve only marginally since the aggre- 
gates do not grow so much. For t > 0.3 Myr the curves overlap 
as the steady-state size distribution has been reached. Despite 
the modest amount of growth, there is still a noticeable increase 
in k of a factor 2-3 for the near- and far-IR opacities. Porous 
aggregates enhance the opacity at these wavelengths. 

In Fig. [6}} the extinction k is plotted for ice-coated distribu- 
tions, corresponding to Fig.QJ). For ice-coated models the size 
distribution evolves significantly (Fig.QJ)), and this is reflected 
in the spectra of Fig. [6J3. Initially, for t < 1 Myr, the trend reflects 
the models without coating, shown in Fig.|6^. The near-IR opac- 
ity significantly increases since aggregates grow to sizes com- 
parable to these wavelengths. Meanwhile, the opacity at wave- 
lengths approaching 0.1 /mi strongly decreases. The opacity is 
especially boosted at the wavelengths that correspond to the 
particle size that contains most of the mass, i.e., that peaks in 
Fig. ID). 

At late times (f > 1 Myr) this trend continues and any spec- 
tral feature is significantly reduced in strength or disappears al- 
together. For the 1 and 3 Myr runs the opacity becomes gray at 
short wavelengths, reflecting the predominance of big, a > 10 
fim, aggregates. However, at t = 10 Myr, there has been some 
replenishment of small grains due to fragmentation among ag- 
gregates (see Fig.QJi), and this is reflected in an increase in the 
opacities at short wavelengths. Also, it is clear that the sub-mm 
opacity has increased significantly due to the existence of ~100 
jum aggregates (Fig.QJi). 

4.2. Implications for r s a vs. E(J - K s ) 



(5) 



In Sect. l3.2l we already discussed the ratio 



q = 



Tsfl 



E(J-K S ) 



0.93- 



Kj ~ KKs 



where t s h is the optical depth in the 9.7 silicate absorption fea- 
ture (i.e., with the continuum subtracted), E(J - K s ) the near- 
IR color excess and ki = r,-/S with S the dust column density. 
This quantity has been used as an indicator of grain growth. 
Observations show that in the outer parts of dense cloud this ra- 
tio is equal to the value observed for the diffuse ISM (q « 0.34; 
iRoche & Aitken1ll984t Iwhittetll2003t fChiar et al]l2607h : but we 
can expect that grain processing (ice mantle formation and ag- 
gregation) will alter the q value. 

Figure [7] shows the temporal evolution of q for the four ag- 
gregate models under consideration. Per model we plot t s u and 
E(J — K s ) for each of the curves shown in Fig. [6] and for each 
of the six times, and connect these points by a line. For the 
(sil , gra) model of Fig. [6^ the points (connected by the black 
line) lie to the right of the ISM-curve. The g-value is lower than 
the ISM, q « 0.20 - a reduction that can be explained as the re- 
sult of type II mixing (see Sect. l3.21 >. In the ice-coated model 
((ic-sil , ic-gra); dashed black line in Fig. [7]) the g-value 
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1CT 1 10° 10 1 10 2 10° 10 1 10 2 10 3 

Wavelength A [jUm] Wavelength A \p.m\ 

Fig. 6. The extinction for the time-dependent porous aggregates distributions at several times. In each model silicates and graphite 
grains are mixed at a 2:1 ratio. See Table[T]and the text for a description. The observed ISM extinction curve (thin solid line) is 
shown for reference. 




5000 10000 15000 

Color excess, E(J — K)fL^ 

Fig. 7. The opacity in the 9.7 fim silicate feature vs. near-IR color 
excess for the four aggregate models presented in Fig. [6] Each 
model is identified by a color/line style as indicated in the leg- 
end. The initial (f = 0.03 Myr) value of r s ;i and E(J - K) per 
unit column density is indicated by an empty square, the final 
value (t — 10 Myr) by an empty circle, and intermediate times 
are indicated by black dots. 



start off from the diffuse ISM line. Initially, the behavior is very 
similar to the sil+gra model: both the E(J - K s ) color excess 
and t s ;i increase and q reaches a minimum value of 0.19 after 
t = 3 x 10 5 yr. However, after t > 3 x 10 5 yr the effect of the 



ice-driven coagulation becomes prominent. At late times (f > 1 
Myr) the line approaches the origin as both indicators vanish. 

Thus, the bare grain (sil.gra) and the ice-coated 
(ic-sil , ic-gra) models succeed in redu cing q below the 
ISM- value, which agrees with the trend of the dChiar et all 2007. 
l201ll) studies. However, it must be recognized that q in these 
data is the net value consisting of two contributions: diffuse ma- 
terial from the cloud at an ISM g-value; and 'processed' material 
from the core. The observed g-value, then, is an upper l imit to the 
g-value of the core material. Indeed. IChiar et al .T (l2007h suggest 
that the q value for the core may be much less than q\sM- They 
speculate that in the core the silicates grains (the carriers of the 
9.7 fim feature) coagulate very rapidly - essential eliminating 
the 10 jjm feature - whereas the carbonaceous grains (the car- 
riers of the near-IR extinction) are not involved in this process. 
Since strong growth is, within the context of our simulations, 
associated with ice-coated grains, this implies that the silicates 
succeed in acquiring ice mantles, whereas the graphite grains 
fail to do so. Although the physical reason for this dichotomy is 
unclear, we next investigate two additional models in which the 
evolution of the silicates and graphites are decoupled. 

First, we consider a scenario where the graphite grains are 
non-coated, but where the sticking behavior is still provided by 
the ice-coated distributions of Fig.QJ). In Fig. [(J the only differ- 
ence from Fig. [(J then is that the graphite contribution present 
in the aggregates is not coated. For this reason the trend in q is 
similar. Note that during the initial phase (f < 3 x 10 5 yr, where 
q » .17) the curve very well follows the trend of Chiar et al. 
(1201 ll) . Like before, this reduction is primarily the result of ice- 
mantle formation and aggregation (e.g., type II mixing), rather 
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than the very growth of these aggregates; there is no strong evo- 
lutionary trend. 

Alternatively, we presume that - again, for unexplained rea- 
sons - silicate and graphite grains fail to interact at all. Then, the 
graphite grains follow the size distribution as given in Fig.QJ, 
whereas the ice-coated silicates follow Fig.QJ). Equivalently, we 
can say that the materials are spatially separated (type I mix- 
ing). As one can see from Fig. [21 this seems to be the only 
way to obtain low q values - provided the timescales are long 
enough (~Myr) to coagulate away the silicate feature. The op- 
tical depth in the silicate feature disappears as the silicates co- 
agulate, whereas the near-IR excess remains large, due to the 
presence of small graphite aggregates. 

4.3. Opacities at sub-mm wavelengths 

Another potential powerful indicator for grain growth is the be- 
havior at far-IR wavelengths. However, in the Rayleigh regime 
(A » particle size) the opacities are not expected to be a strong 
function of the size of the particles. Therefore, we see initially 
(f < 1 Myr) only little variation among the models; and this vari- 
ation can best be attributed to the way the components are mixed 
and whether ice-coating is involved. Table [2]provides the opaci- 
ties at 150, 350, and 500 /mi - the wavelengths corresponding to 
the SPIRE bands - and at 850 /mi for the four aggregate models 
we considered here. These opacities are given for several coag- 
ulation timescales. 

The differences in opacity between the four aggregate mod- 
els after 10 5 yr reflect the uncertainty in the initial conditions of 
the dust distribution (i.e., ice-coated or not) rather than indicat- 
ing grain growth. After t — 10 6 yr very little has changed. One 
interesting point is the different behavior seen in Fig. [6}; and d, 
i.e., between the (ic-sil,gra) and ic-sil+gra models. In 
the former the sub-mm opacities decrease by about a factor two, 
whereas in the latter there is an increase by a factor two (see also 
Fig ©- This difference can be attributed to the applied mixing 
rule: the IR-opacity of free-floating graphite aggregates strongly 
increases with (modest) aggregation; however, if graphite is em- 
bedded in aggregates that are dominated by ice-coated silicates, 
its signature is suppressed. 

However, after 10 7 yr of coagulation there is a strong in- 
crease, up to a factor 20, in the opacity of the ice-coated silicate 
models, which is the result of grain growth. Most of the mass 
is now in 100 /mi size particles (see Fig.QJi) and this boosts the 
emissivity at sub-mm wavelengths. Although a factor of 5 dif- 
ference in k may be explained by the model setup or material 
properties, only strong growth up to 100 /mi sizes can explain 
opacity enhancements by more than a factor of 10. 

In Table [2] we ha ve compared our op a city c alcula- 
tion with those of Ossenkopf & Henning (1994). The 



lOss 
rfl99' 



Ossenkopf & Henning (1994) are based on an aggregation 
scenario that features only 'hit-and-stick' coagu lation, meanin g 
that large fluffy (fractal) structures will form (Ossenkopf 19931) . 
This, as well as many other factors, are different from our setup. 
Nonetheless, after 10 5 yr of coagulation, the sub-mm opacities 
are comparable, which gives support to both approaches. 
However, the high opacities after t — 10 7 y r are not present 
in any of the Ossenkopf & Henning (1994) tables, also not 
at the higher gas densities n g they consider^] Generally, one 
can expect that a higher gas density (which means a shorter 



collision timescale) should r esult in faster growth, but in the 
lOssenkopf & Henningl |[l994) calculations, this trend is not so 
obvious (see also the discussion in Paper I). 

An increase in the sub-mm slope has seen recent interest 
with Pl anck observations of cold galacti c clouds indicative of 
P > 2 dPlanck Collaboration et alJlioi lalTbh. Likewise, detailed 



1 Note that the most f requently quoted opacities from the 
lOssen kopf & Henning] i 19941) study correspond to the fifth column of 
the opacity table in their paper (OH5), as calculated from a gas density 



studies of in dividual cores (e.g. JSchnee et alj2010tlShirlev et all 
I2005L 1201 ll) indicate that /3 is likely to be larger than 2. A fa- 
vored explanation for the high /^-values in these cold clumps is 
that the opacity becomes temperature dependent. Several labora- 
tory studies show that lo w temperatures suppress the emissivity 
at longer wavelengths dMennella et alJll998l:lBoudet et~ail2 005: 
Demyk et al., in pr ep). The observed c orrelation between /3 and 
decreasing T (e.g., Des ert et al.l 12008) gives credibility to this 
scenario. Dust aggregation is not a requirement. 

On the oth er hand, aggregation can also change the ft value 
dWrightHl987l) . In our study, the formation of ~100 /mi aggre- 
gates boosts the emissivity at the short end of the sub-mm range 
(~100 //m; see Fig. [TJ, thereby increasing ft. Thus, the high ft re- 
sults from an increase in the sub-mm emissivity - not a decrease, 
as in the low T scenario. We may therefore distinguish the two 
scenarios by comparing the derived sub-mm opacity to the opac- 
ities at other wavelengths. In the aggregation scenario a larger ft 
also results in a larger Ksub-mm^near-iR ratio (see Table|2]i. We do 
not expect this trend to emerge, however, when the increase in ft 
is driven by a low temperature. In this light, it is interesting that 
our findings regarding the opacity ratio between 850 and 2.2 /mi 
for the (ic-sil,ic-gra ) model after 3 x 10 6 yr are consistent with 
IShirlev et ail (120111) for the B335 class source. 

Polarization studies, finally, can also provide valuable in- 
sights about the typical size of aggregates. Small aggregates 
can be highly anisotropic, even while the aggregation pro- 
cess is isotropic, whic h results in their alignment in a field 
dBotet & Rannouir2003l) . However, with increasing size the op- 
tical anisotropy is expected to diminish. 



5. Summary and Discussion 

In this study we have investigated the effects of grain growth 
on the mass-weighed opacities for dense molecular clouds, fo- 
cusing on IR w avelengths. From previous works (Paper I and 
iMin et al]|2008l) we have adopted a time-dependent size distri- 
bution and a method to calculate the opacities using an effective 
medium approach. In our analysis we have mainly focused on 
indicators for grain growth: the ratio between the IR-colors to 
the strength of the 9.7 //m silicate feature and the behavior at 
sub-mm wavelengths. These we have attempted to link to obser- 
vational studies. 

In all of this, it should be emphasized that the opacity values 
present in this work flow from a long chain of modeling, during 
which several assumptions have been applied. For example: the 
collisional model and the corresponding aggregate size distribu- 
tions of Paper I are based on the approximation of similar-size 
grains (0.1 //m in this case); the gas density has been fixed at 
«gas = 10 5 cm~ 3 ; and we only considered silicate, graphite and 
ices as material properties. Varying these parameters will give a 
different quantitative picture. The merit of this work, then, pri- 
marily lies in the trends that are exposed when ice-mantle and 
aggregate formation/growth are taken into account. 

Still, the current setup already provides a considerable 
amount of freedom to tune the opacity values. Figure [4] shows 



of 10 s g cm 3 . The n g = 10 5 data in Table [2]has been retrieved from the 



online data. 
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Name 


Time 




Opacity [cm 2 g '] 




^850 '■ K2.2 


P 






150 yum 


350 iim 


500 pm 


850 yt/m 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(sil,gra) 


TrX 


49 


8.6 


4.4 


1.7 


3.1 x 10~ 4 


1.8 




10 6 


56 


9.6 


4.8 


1.7 


2.5 x 10~ 4 


1.9 




3 x 10 6 


56 


9.5 


4.8 


1.7 


2.5 x 10~ 4 


1.9 




10' 


56 


9.5 


4.8 


1.7 


2.5 x 10~ 4 


1.9 


(ic-sil,ic-gra) 


10- 


50 


6.8 


3.2 


1.1 


3.0 x 10~ 4 


2.0 




10 6 


41 


5.5 


2.7 


0.92 


7.9 x 10 5 


2.0 




3 x 10 6 


120 


10 


4.0 


1.2 


4.6 x 10~ 4 


2.3 




10' 


410 


150 


58 


10 


7.2 x 10~ 3 


3.3 


(ic-sil,gra) 


10 


60 


10 


5.1 


i .y 


3.0 X 10 


1.8 




10 6 


36 


5.2 


2.5 


0.87 


7.7 X 10~ 5 


2.0 




3 x 10 6 


120 


10 


3.9 


1.1 


4.8 x 10 


2.4 




10' 


370 


170 


69 


12 


7.3 x 10~ 3 


3.4 


ic-sil+gra 


10 5 


42 


5.3 


2.5 


0.86 


1.2 x 10~ 4 


2.0 




10 6 


64 


9.1 


4.6 


1.8 


1.3 x 10~ 4 


1.8 




3 x 10 6 


100 


11 


5.3 


1.9 


1.8 x 10~ 4 


1.9 




10 7 


290 


75 


29 


5.9 


6.3 x 10- 4 


3.0 


OH94, no ice 


10 1 


24 


6.5 


3.9 


1.6 


2.5 x 10~ 4 


1.6 


OH94, thin ice 


10 5 


35 


7.8 


3.9 


1.4 


7.2 x 10- 3 


1.9 


OH94 thick ice 


10 5 


57 


10 


5.0 


1.8 


7.3 x 10" 3 


1.9 



Note. — Opaciti es at sub-mm wav elengths for the four aggregate model at several times (Cols. 1-6). The bottom rows give the corresponding 
opacities from the Ossenkopf & Henning ( 1994) models for a gas density of n glis = 10 5 cm 4 . (Col. 7) Ratio between the 850 yum and the 2.2 //m 
opacity. (Col. 8) Index of the sub-mm power law slope, k a cc with /? derived from the opacity values at 850 and 500 /jm. 



that a considerably variation in opacity can already be achieved 
by merely varying the way to mix silicate and graphite grains 
(or their ice-coated variants), without taking into account the ef- 
fects of (strong) aggregation. These variations affect the ratio 
between the 9.7 fim silicate feature and the near-IR color excess, 
q = T s ji /E(J-K), which has been found to be constant in diffuse 
clouds, q ss 0.34, but deviating to lower values for dense cores. 
Then, if we fix the ratio silicates: graphite such that q « 0.34 for 
the diffuse ISM, we have found that an ice coating increases q if 
the grains are spatially mixed. However, if the materials involved 
are mixed on the level of the aggregate the g-value substantially 
decreases below the diffuse ISM line. 

The effects of coagulation add to these findings. Its effect 
can best be understood in terms of the size of the most domi- 
nant dust aggregates in the distribution that evolve as function of 
time. First, the near-IR color excess increases as aggregates grow 
to ~yizm sizes. Further growth requires grains to be ice-coated. 
The near-IR color excess then first decreases as the dominant 
aggregate size becomes larger than the size corresponding to the 
J and K wavelength bands. The strength of the 9.7 fim feature 
decreases only thereafter when most of the dust mass becomes 
locked up into larger and mo re compact aggregates. Thus, the 
finding of Ivan Breemen et al.l (1201 lh that for their sources it is 
the J - K excess that affects q, implies, within the context of this 
work, that their objects are only moderately affected by coagu- 
lation. 

The same conclusion may be drawn from the observed ra- 
tio of the strength in the 9.7 um feature vs. near-IR color ex- 
cess, as presented by IChiar et al.l d2007ll20Tlh . This ratio (q) is 
lower than the ISM value by roughly a factor of two, and this 
is in good agreement with our ice-coated models for a timescale 
t < 3 x 10 5 yr. We find that this decrease is, however, primarily a 
result of the formation of aggregates that consist of different ma- 
terials, rather than their growth. Indeed, for stronger growth, we 
would expect larger q values to reemerge. To obtain even lower 
^-values - applicable to, possibly, the dense cores - requires a 



decoupling of the carbonaceous and silicates, in the sense that 
the former species does not accrete ice mantles, whereas the lat- 
ter does, and where interaction between the carbonaceous and 
silicate species is inhibited. Under these conditions, our simula- 
tions show that q <sc 1 is possible, since the carbonaceous grains 
only mildly coagulate, whereas the silicates aggregates become 
large. Although such a scenario may seem a bit far-fetched, it 
is the only way to produce q <K 1. It is therefore very helpful 
to further ascertain observationally the state of the dust in these 
high density cores. 

In order to better constrain the evolutionary state of a cloud, 
multi-wavelengths observations are essential, e.g., from the 
shape of the silicate feature, the scattering behavior at mid-IR 
wavelengths, or from sub-mm data. As a general rule of thumb 
we have found that the opacity peaks for the wavelengths cor- 
responding to the size of the dominant aggregates in the distri- 
bution. This behavior is perhaps most obvious at sub-mm wave- 
lengths. Initially, little variation is expected (whereas the near- 
IR color excess, for example, are already strongly affected). 
However, if grain growth yields aggregates of size ~100 fim, the 
sub-mm opacities are strongly boosted and the sub-mm power- 
law index f3 will increase (i.e., the opacity curve steepens). 

Acknowledgements. The authors appreciate the helpful comments of the referee 
and editor. 
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